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We calculate the QED corrections to deep inelastic scatter- 
ing with tagged photons at HERA in the leading logarithmic 
approximation. Due to the special experimental setup, two 
large scales appear in the calculation that lead to two large 
logarithms of comparable size. The relation of our formalism 
to the conventional structure function formalism is outlined. 
We present some numerical results and compare with previous 
calculations. 

PACS number(s): 13.60.-r, 13.60. Hb 



I. INTRODUCTION 

One of the major tasks of the experiments at the 
HERA collider is the determination of the structure func- 
tions of the proton, F 2 {x,Q 2 ) and Fr / (x,Q 2 ), over a 
broad range of the kinematical variables. The exten- 
sion of these measurements to the range of small Bjorkcn 
x < 10 -4 and Q 2 of a few GeV 2 is of particular inter- 
est, because it will help improve our understanding of the 
details of the dynamics of quarks and gluons inside the 
nucleon 0. 

In order to separate Fl(x,Q 2 ) from F2(x, Q 2 ), it is 
necessary to measure the cross section for ep — » eX with 
different center-of-mass energies. However, instead of 
running the collider at reduced beam energies, one can 
employ a method suggested by Krasny et al. that uti- 
lizes radiative events. This method takes advantage of a 
photon detector (PD) in the very forward direction, as 
seen from the incoming electron beam. Such a device 
is part of the luminosity monitors of the HI and ZEUS 
experiments. 

The idea of the method is that emission of photons in 
a direction close to the incoming electron can be inter- 
preted as a reduction of the effective beam energy. The 
effective beam energy for each radiative event is deter- 
mined from the energy of the hard photon that is ob- 
served in the PD. In fact, radiative events were already 
used in a measurement of the structure function F 2 down 
to Q 2 > 1.5 GeV 2 §. 

The possibility to use radiative events for structure 
function analyses was already discussed by Jadach et 
al. H. However, these authors used radiative events with 



tagged photons in order to reduce the radiative correc- 
tions and to simplify the F 2 analysis. They did not con- 
sider higher order corrections explicitly, and it is not clear 
whether their method allows an extraction of Ft, with- 
out running at lower beam energies. On the other hand, 
the feasibility of a quantitative measurement of Ft, at 
HERA using the method of @, for Q 2 below 5 GeV 2 and 
x around 10~ 4 , has been studied in || and considered 
possible. 

A more general treatment of the Born cross section 
for ep — > ejX that is also valid for non-collinear photon 
emission and in the range of high Q 2 can be found in the 
paper by Bardin et al. S . 

For an accurate description of the corresponding cross 
section one has to consider the radiative corrections. In 
this letter we calculate the QED radiative corrections to 
deep inelastic scattering with one tagged photon to lead- 
ing logarithmic accuracy using the structure function for- 
malism . We will in particular discuss the appearance 
of two large logarithms corresponding to two large scales, 
which are due to the experimental setup peculiar to the 
HERA experiments. Finally, we will show some numeri- 
cal results and compare our findings with the calculation 
by Bardin et al. M. 



II. KINEMATICS AND LOWEST ORDER CROSS 
SECTION 

To begin with, let us start with a brief review of the 
kinematics adapted to the case of deep inelastic scatter- 
ing with one exclusive hard photon, 



e(p e ) + p(P) - e(p' e ) + 7 (fc) + X(P'), 



(1) 



where the polar angle $ 7 of the photon (measured with 
respect to the incident electron beam) is assumed to be 
very small, i? 7 < $o, with $o being about 0.45 mrad in 
the case of HI. 

A convenient set of invariants that takes into account 
the energy loss from the collinearly radiated photon is 
given by [§: 

Q 2 = -(Pe-p' e ~k) 2 , 
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2P • (p e - P ' e - ky 
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P-(p e -k) ■ 



(2) 



Since we restrict ourselves to collinear photons, it is suit- 
able to parameterize the energy of the radiated photon, 
E 7 , with the help of 



E e — E 1 
EL 



xyS' 



with S = 2p e ■ P. (3) 



The differential cross section for the process ep — > ejX, 
integrated over the photon emission angle < i? 7 < $0j 
reads [01 



a fBorn 

dxdQ 2 dz 
where 

cr (x,Q 2 ;z) 
with 

y 

and 
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with 



i = lnCo, Co 



E 2 Jl 



(4) 

F 2 (x,Q 2 ), (5) 

(6) 
(7) 

(8) 



Here we have neglected terms of order G($q) as well as 
terms of order O^q 1 ), and we have neglected the contri- 
butions from Z° exchange since we are mainly interested 
in the region of small Q 2 (see also ||). R is defined as 
the ratio of cross sections for longitudinal to transverse 
photons, 



III. RADIATIVE CORRECTIONS 

There are two large sources of contributions to the 
radiative corrections from higher orders that contribute 
terms of the order of a/ir ■ hiQ 2 /m 2 , with mj being the 
mass of a light fermion. One type of contributions are 
the corrections to the propagator of the exchanged boson 
between the electron and the hadronic system, while the 
others are due to radiation of photons off the incoming 
electron line, which we will treat in the structure function 
formalism. 

We shall assume a calorimetric experimental setup, 
where a hard photon radiated collinearly to the outgoing 
electron line cannot be distinguished from a bare out- 
going electron, so that final state radiation can be ne- 
glected to the desired leading logarithmic accuracy in ac- 
cordance with the Kinoshita-Lee-Nauenberg theorem || . 
We also assume some minimum experimental cut on the 
transverse momentum of the outgoing hadronic system, 
in order to suppress the contribution from QED Compton 
events [To| ] for a leptonic measurement of the kinematic 
variables. 



A. Vacuum polarization 

The first important type of contribution to the ra- 
diative corrections comes from the vacuum polarization, 
which amounts to replacing the coupling constant a in 
the hard cross section Co (eq. by the QED running cou- 
pling a(—Q 2 ). For the contribution from lepton loops we 
use the well-known one- loop perturbative result. Since 
we are interested in the region of rather low Q 2 , we use a 
parameterization by Burkhardt and Pietrzyk [ pd| for the 
hadronic contribution. 



R{x,Q 2 ) 



F L 



(9) 



The allowed range for z follows from (g) and the restric- 
tion < y < 1, 



xS 



< z < 1. 



(10) 



Note that the entire ^-dependence of Co occurs only via 

Eq. (HI) agrees with the Born cross section given in 
H when restricted to the kinematical situation under 
consideration.^] 



1 We used the same set of kinematical variables as in j2| and 
[H . For details on the relation between these two sets we refer 
the reader to the discussion in 101. 



B. Photonic corrections 

In order to illustrate the calculation of radiative cor- 
rections to the cross section ([|) due to photon emission in 
the framework of the structure function formalism HQ , 
let us first note that this cross section is already a QED 
correction: it is that part of the inclusive first order cor- 
rection to deep inelastic scattering (ep — > eX), which 
is selected by requiring an exclusive hard photon seen 
in the PD. This may be exhibited by writing down the 
radiatively corrected inclusive cross section to DIS as a 
convolution of the electron non-singlet structure function 
-Dns(z, Q 2 ) with the hard cross section (||) 



d o-Rc 
dxdQ 2 



dz Dm(z,Q 2 )a (x,Q 2 ;z). 



(11) 



The electron non-singlet structure function, 



2 



D m (z,Q 2 )=S(l-z) + ^LP^(z) 

Z7T 



(£) 2 ^ 2 Hz) + 0((aLf), (12) 



depends on the large scale Q 2 only via the large logarithm 
L = InQ 2 /to 2 . It is known to properly sum the leading 
contributions (aL) n to all orders in perturbation theory 

0- 

We now give the relevant coefficients of the power series 
expansion of -Dns ■ Introducing a small auxiliary parame- 
ter A that serves as an infrared (IR) regulator to separate 
virtual+soft and hard photon contributions, the first two 
coefficients of the expansion of -Dns are 

P^' 2 \z) = P 5 (1 ' 2) • (5(1 - z) + Pg' 2 \z) • 6(1 - A - z), 
with 



1 + z A 
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(13) 



and similarly (see e.g. M), 
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(14) 



Inspecting (0), it is obvious that the logarithmic piece 
of P(z) is contained in the first order correction to the 
inclusive cross section; the difference in the logarithms 
(L — Lq) is accounted for by the remaining phase space 
due to photons emitted with angle larger than z9 - 

Let us now turn to the contributions from the second 
order. Since the maximum emission angle i9 is about 
0.45 mrad for the HERA detectors, the "large logarithm" 
Lq appearing in (Q) turns out to be moderate, 



Co « 600, L « 6.4 

so that the complement 

L\ = L — Lq w 6.5 . 
(for e.g., Q 2 = 0.1 



at E e = 27. 5 GeV, 



15.7 

.1000 GeV 2 ) 



is of similar magnitude or even larger than Lq. For this 
reason we have two large logarithms Lq, L\ entering the 
game. 

We shall separate the contributions into vir- 
tual+collinear, soft+collinear, two collinear photons, and 
one collinear plus one semicollinear photon. 

The sum of the contributions of virtual+soft correction 
to collinear photon emission, with the emission angle of 
the soft photon being integrated over the full solid angle, 
but the hard photon only over the angular range of the 
PD, can be obtained from the expression for the one- loop 
Compton tensor. It reads |L2] 



Lq 



(Lq + L 1 )P( 1) - (Lq 

2 ;z). 



2Li)lnz 



(15) 



The contribution from two hard photons in the PD, 
with the energy fraction of each being larger than A, is 
found to be 13 



/ a 
V2^ 



2 M 

2 

xa (x, 



P.. 



2P, 



P> 



(i) 



In z 



(16) 



with 2 = 1— (^2 E^/Ef., where E 7 is the total photon 
energy registered in the PD. 

Finally, we consider the contribution from one collinear 
photon that hits the PD, while the other one is emitted 
at an angle larger than t?o- 

Let us investigate the regions of phase space where the 
large logarithmic contributions originate. For the hard 
photon that hits the PD, the major contribution comes 
from the region of polar angles ~ m e /E e <g; ^o- 
Similarly, for the other photon the biggest contribu- 
tion comes from polar angles close to the lower limit, 

(2) 

^7 ^ $o- Thus, the leading contributions come essen- 
tially from the region where (tf^M^) 2 w Co > 1. One 
can easily see that one obtains the double logarithms en- 
tirely from the contribution where the photon hitting the 
PD is emitted first, with the "lost" photon being emit- 
ted second, while the reversed case is suppressed if the 
following condition is satisfied: 



X2 ■ Co > 1) 



E {2) 

where X2 = 

E K 



(17) 



Thus we shall assume in the following A <C 1, but A -Co ~S> 
1. 

Taking this ordering of photon emissions into account, 
the contribution from one tagged plus one undetected 
photon can be calculated by means of the quasireal elec- 
tron method [G3 



/ a 



v .,_, L Q L l -P^(z) 



A z V * 



xa (x,Q 2 ;z - x 2 ), 



(18) 



where &q is understood to be expressed by the kinematic 
variables x, Q 2 , of the hard subprocess, 

a (x, Q 2 ; z - x 2 ) = <7q(x, Q 2 , z - x 2 ) • J(x, Q 2 ;x 2 )- (19) 

Note that this contribution explicitly depends on the ex- 
perimental determination of the kinematical variables x 
and Q 2 , since the almost collinear emission of the sec- 
ond photon shifts the "true" kinematical variables (£, Q 2 ) 



3 



with respect to the measured ones (z,Q 2 ).g The Jaco- 
bian J accounts for scaling properties of the chosen kinc- 
matical variables under radiation of the second photon. 
The upper limit of the ^-integration in (|l8| ) is given by 
either some experimental cut on the maximal energy of 
the second photon, or by the kinematical limit, which also 
depends on the choice of the experimental determination 
of the kinematical variables, as we will discuss later. 

After change of variables x 2 = zu, uq = a;™ ax /z, the 
integral in ( |l8| ) may be conveniently decomposed into 
IR divergent (A dependent) and IR convergent contribu- 
tions as (suppressing the arguments x, Q 2 ) 



1 + ^ho = 

aLi 

1 



a(-Q 2 ) 
a(0) 



{+ aLo lz± P m {z) 



+ 



duP£\l~u) 



4tt 1 + z 2 6 
a (z(l - u)) 



<tq(z) 



- 1 



[ duP«(l-u) \ 

J Uq J J 



(22) 



Note that the contribution from the undetected hard pho- 
ton depends on the choice of kinematical variables and 
on the upper limit uq. 



A 



dx 2 D (l) (, X A~( \ 

Pq Ml a (z - x 2 ) = 



du Pg\l - u) [(a (z(l - u)) - a {z)) + a {z)} 



A/z 

ctq(z) 



duP£\l-u 







a (z(l - u)) 



+2\nz-P, 



(i) 



/ duPg } (l-u) 

Ju 



, (20) 



where in the last step we have extended the u-integration 
of the IR convergent piece to due to the smallness of 
A, and we have used the property 



du P (1 \u) = 



in the simplification of the IR divergent piece. 

If we sum the contributions (|l5|), (|l6|), and (p"8|), we 
see that the dependence on the auxiliary parameter A 
cancels, as it should. 



IV. RESULTS FOR HERA 

In the presence of the an undetected (lost) photon, the 
determination of the kinematical variables x, Q 2 becomes 
ambiguous, it will de pen d on the method chosen, and in 
turn the corrections ( P2[ ) will depend on this choice. 

For HERA, several methods are being used to deter- 
mine the "true" kinematical variables x,Q 2 , in order to 
reduce systematic errors or to control the influence of 
initial state radiation (see e.g., [|l5| for an illuminating 
discussion of the latter). For the sake of brevity we will 
discuss only the so-called electron method (E) and the 
Jacquet-Blondel method (JB). 

In the case of the electron method, the kinematical 
variables (0) are determined via 



Qi = 4:El n E' e cos^ e /2, 
E' 

He = 1 - -r^s sin 2 i3 e /2, x e 



with 



Ef - y e s^' 



s off = 4 E f E E f =E E zE 



(23) 



C. The radiatively corrected cross section 

We may now write down our final result for the ra- 
diative corrections. Since we restricted ourselves to the 
leading logarithmic corrections, it follows that the radia- 
tively corrected cross section may be written in factorized 
form, 



d ctrc 



d d cr B o 



da;dQ 2 dz dxdQ 2 dz 



{l + 5 ho {x,Q 2 lZ )) . (21) 



where we retain in the correction factor (l + ^ho) only the 
logarithmic terms Lq,Li from ([To]), jlq), and (|l8|), and 
the contribution from vacuum polarization, 



and E 1 is the energy deposited in the PD. Radiating an 

(2) 

additional almost collinear photon with energy Ey = 
x 2 E e leads to a shift of the "true" variables x,Q 2 , of 
the hard subprocess with respect to the measured ones, 



Q z = Ql- z', y 



Ve+Z'-I 



Ve + Z 



ZT — i"' ( 24 ) 
z' — 1 



with z' = 1 — x 2 /z = l — w. The kinematical limit for 
the undetected photon follows from x < 1, 



z' > 



1 ~ Ve 

1 - x e y e 



(25) 



The Jacobian (^9|) for this choice of kinematical variables 
is 



2 For a recent review on the influence of initial state radi- 
ation on the experimental determination of the kinematical 
variables see e.g., Cq|, and references cited therein. 



J(x e ,Ql)=( VeZ 

\y e + z' - 1 



(26) 



In the case of the Jacquet-Blondel method, the kine- 
matical variables are determined using 
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2/JB = 



<3jB 



2Ef ^ 

Erf 
i - 2ajb ' 



[E h -p z ,h), 



Q 



JB 



yjBS 



off : 



(27) 



where Eh,p z ,h,Ph are the energy, longitudinal and trans- 
verse momentum components of the final state hadrons. 

The relation between the measured and "true" vari- 
ables in presence of an undetected photon are now given 
by 



-.2 _ 1 ~ 2/JB 



z> • 



f 



'1 / / ' ; 

1 - ym/z' 

with the Jacobian 

J{*jb,C%b)=( , 1 ~ VJB / , 
V 1 - ym/z' 

and with the kinematical limit being 

yjB 



x.m- 



ym 



ym/z' 
(28) 



(29) 



z' > 



1 - x m {l - y JB ) ' 



(30) 



Comparing (24) and (|28|), one finds that one can in 
principle determine the energy of the undetected photon 
(assuming that it is emitted almost collinearly) via 



z = 1 + Vjb - ye- 



(31) 



This may be used to impose an experimental cut to re- 
duce the size of the radiative corrections. 

Let us now present some numerical results for the ra- 
diative corrections (^2|). As input we used 

E e = 27.5 GeV, E p = 820 GeV, t9 = 0.45 mrad, 

(32) 

the GRV94-LO structure functions Q from PDFLIB Q 
for x > 10" 3 , the BK parameterization^ @ for GRV94 
for x < 10~ 3 , and for simplicity a fixed value R — 0.3 (see 
e.g., || for a table of measured R values). Figure 1 shows 
the corrections 8\ w for an energy of 5 GeV deposited in 
the PD. 



3 We have verified that the results for the corrections do not 
change significantly when we use the ALLM [JOI parameteri- 
zation for x <y 10 -3 , although the difference in the predictions 
for the low-a;, low-Q 2 range affects the corrections already for 
x around 1CP 4 . 
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FIG. 1. Radiative corrections <5ho (|2^) in leptonic variables 
for different values of x and a tagged photon energy of 5 GeV. 

At small y e , the corrections are negative, since the con- 
tribution from virtual+soft corrections dominates, be- 
cause the kinematical limit on the energy of the unde- 
tected photon tends to zero as y e — > 0, 

y e (l-x e ) 



"0 = 1 



(33) 



nn 1 - x eVe 
For large y ei the phase space for photon emission in- 
creases, increasing also the shifts between "true" and 
measured variables (g4j), which leads to large positive 
corrections as y e — > 1. 
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FIG. 2. Radiative corrections <5i 10 ( p2| ) in Jacquet-Blondel 
variables for different values of x and a tagged photon energy 
of 5 GeV. 
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Figure 2 shows the corrections <5ho for Jacquet-Blondel 
variables, again for a tagged photon energy of 5 GeV. 
In this case the corrections are negative for j/jb — > 1, 
because in this limit the phase space for the undetected 
photon tends to zero, 



Uq = 1 



(1 - 3Mb)(1 - 2/jb) 
1 - x JB (l - j/jb) 



(34) 



whereas for small ?/jb the corrections remain moder- 
ate. Since the Jacquet-Blondel variables correspond to 
a "more inclusive" measurement than the leptonic vari- 
ables, the corrections due to radiation of an additional 
photon are generally smaller. 



in ||. More details on this interference between photon 
emissions can be found in |jl3|| . Furthermore, it seems 
that the authors omitted the statistical factor 1/2!, while 
summing the contributions (iib, iic) of the semicollincar 
kinematics. 

Some care is needed in the comparison of numerical 
results since Bardin et al. choose kinematical variables 
different from ours. Taking into account the abovemen- 
tioned factor of 2 and the difference due to our treatment 
of the ordering condition for the photon emission, we can 
qualitatively reproduce their results as presented in fig- 
ure 6 of their paper in the range of smaller y values (thus 
confirming the lack of the symmetry factor). 



V. DISCUSSION AND CONCLUSIONS 

In this letter, we have studied the radiative correc- 
tions to deep inelastic scattering with tagged photons at 
HERA taking into account only the leading logarithms. 
The relevant expansion parameter in the present case is 
not simply a/ir ■ L, but we have two large logarithms 
Lo, L\ = L — Lo which happen to be of similar order of 
magnitude due to the particular choice of the geometri- 
cal acceptance of the photon detector. The corrections 
do depend on the choice of the experimental determina- 
tion of the kinematical variables, but they turn out to be 
of the order of 20-40% in most regions of phase space for 
leptonic variables, and below 10% for Jacquet-Blondel 
variables. Large negative corrections occur in those re- 
gions of phase space that are dominated by soft photon 
emission. These large corrections are however well under 
control once one takes into account the resummation of 
multiple photon emission |||7]]. 

We have not considered the contributions from next- 
to-leading order, (a/ir) 2 x L, which may be sizable since 
the individual logarithms Lq, Li are not very large, espe- 
cially for the experimentally interesting region of low Q 2 . 
However, the calculation of those terms is more involved 
and will be presented elsewhere [ pp| . 

We have compared our results with the results given 
in ||. At the Born level, we find very good numerical 
agreement between (^J) and the program HECTOR |2l|] if 
Z-exchange is neglected, which is a good approximation 
since the experimental analysis is restricted to the region 
of not too large Q 2 ||. However, at the level of radia- 
tive corrections, there are major differences. First, we 
note that our derivation of the calculation of the correc- 
tions that is based on the structure function formalism 
|j|,[7| disagrees with formula (4.1) given in their paper. 
In particular, the coefficient of the leading logarithmic 
term derived from (4.1) equals two times our coefficient. 
The naive usage of the quasireal electron approximation 
in their formula leads to a loss of interference of ampli- 
tudes, as was claimed by the authors. However, the inter- 
ference actually does contribute leading logarithms when 
both photons hit the PD, contrary to the statement given 
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